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Reduction of boundary effects in the spiral MRI experiment PROMISE 
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Magnetorotational instability (MRI) is one of the most important and most common instabilities in astrophysics. It is 
widely accepted that it serves as a source of turbulent viscosity in accretion disks - the most energy efficient objects in the 
Universe. However it is very difficult to bring this process down on earth and model it in a laboratory experiment. Several 
different approaches have been proposed, one of the most recent is PROMISE (Potsdam-ROssendorf Magnetorotational 
Instability Experiment). It consists of a flow of a liquid metal between two rotating cylinders under applied current-free 
spiral magnetic field. The cylinders must be covered with plates which introduce additional end-effects which alter the 
flow and make it more difficult to clearly distinguish between MRI stable and unstable state. In this paper we propose 
simple and inexpensive improvement to the PROMISE experiment which would reduce those undesirable effects. 
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1 Introduction 

Velikhov (1959) showed that, for ideal magnetohydrody- 
namics, an axial magnetic field applied to a flow of a liquid 
metal between two concentric, differentially rotating cylin- 
ders (Taylor-Couette flow) can lower the critical rotation 
ratio or even destabilize the flow, although it is hydrody- 
namically stable (when the Prandtl number is large enough). 
That type of instability is called magnetorotational instabil- 
ity, and not long ago Balbus & Hawley (1991) demonstrated 
that it plays an important role in astrophysics. MRI serves as 
an essential mechanism for transporting angular momentum 
in a wide range of astrophysical objects, stellar interiors, jets 
and in particular it is crucial for process of accretion where 
it provides necessary amount of turbulent viscosity. Recent 
experiments by Ji et al. (2006) suggested that purely hydro- 
dynamical quasi-Keplerian flows, i.e. Taylor-Couette flows 
which resemble Keplerian disks, do not provide viscosity 
required to transport the angular momentum effectively, and 
therefore this instability is ruled out as a source of viscous 
turbulence in the disks. 

One of the most convenient laboratory models for MRI 
is still the magnetohydrodynamical, cylindrical Taylor- 
Couette flow with an imposed external magnetic field along 
the axis (Riidiger & Zhang 2001 ; Ji, Goodman & Kageyama 
2001). For magnetic Prandtl numbers Pm w 1-10 the ax- 
ial field (when not too strong) reduces the critical Reynolds 
number for instability and, more importantly, also intro- 
duces instability for flows which are always stable for non- 
magnetic cases (see e.g. Riidiger, Schultz & Shalybkov 
2003; Barenghi et al. 2004). The unstable state is charac- 
terized by classical Taylor vortices. Regrettably, laboratory 
liquid metals possess very small magnetic Prandtl numbers 
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(due to low conductivity) so that for a purely axial field the 
critical Reynolds number is of order O(10 6 ), and therefore 
vast rotation rates are necessary for MRI to grow. 

Recently it has been shown by Hollerbach & Riidiger 
(2005) and Riidiger et al. (2005) that a current-free external 
azimuthal magnetic field in addition to the axial one can re- 
duce the critical Reynolds number to O(10 3 ) which makes 
it much easier to design an MRI experiment. Moreover due 
to symmetry-breaking there exists a drift of Taylor vortices 
associated with the configuration of the applied magnetic 
field. The frequency of the traveling wave that can be mea- 
sured and compared with theory is an important feature of 
this type of instability. 

The idea of an additional toroidal field was success- 
fully implemented in the PROMISE experiment by Stefani 
et al. (2006) where modes corresponding to so called spi- 
ral (or helical) MRI were observed for the first time (see 
also Riidiger et al. 2006; Stefani et al. 2007). Results of this 
experiment also show that in the basic stable state, without 
any toroidal field, there exists a nonzero axial velocity field 
which arises due to presence of the rigid endplates enclos- 
ing the cylinders. These plates, undoubtedly present in any 
real experiment, are responsible for additional effects which 
do not take place for an idealized infinitely long container. 

The boundary layer which exists in the vicinity of the 
endplates consists of an Ekman layer which is the result 
of the rotation of a rigid surface, and a Hartmann layer 
which develops when a conducting fluid is used and an ex- 
ternal axial magnetic field is applied (see e.g. Ekman 1905; 
Roberts 1967). Consequently, the global properties of the 
flow change when compared to the idealized case of in- 
finitely long cylinders: a secondary flow, i.e. two large Ek- 
man vortices appear, and the Hartmann current is drawn into 
the bulk of the fluid. All these effects depend on the mechan- 
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ical and magnetic properties of the lids. In the PROMISE 
experiment one of the lids is made of copper and is attached 
to the outer cylinder, the other one is a stationary plexiglass 
plate. 

In this work we review simple improvements which can 
reduce undesired effects induced by the lids and provide 
therefore the possibility to distinguish more clearly between 
stable and unstable states of MRI. 

2 Numerical model 

We consider two concentric cylinders with radii Ri n = 4 
cm, i? ou t = 8 cm and height H = 40 cm which ro- 
tate with angular velocities Vt m: f2 ut- The rotation ra- 
tio \x = fi out /£!i n is chosen in such a way that the flow 
is hydrodynamically stable, i.e. Rayleigh stability criterion 
d r (r 2 fl) 2 > 0, r being the distance from the axis of rota- 
tion, is fulfilled. Through this paper we use jl — 0.27. For 
infinite cylinders the basic rotational profile for the flow is 
the Couette solution 

n Q (r)=a+^, (1) 

where a, b are constants dependent on radii and rotation 
speeds 



a = fii, 



fj 2 



(2) 



1 - f] 2 - 1 1 
The external magnetic field (steady, current free) has the 



form of 
Bo = B a 

and its strength is measured by the Hartmann number 



(3) 



Ha = B Q J Rini - Rout Rin \ (4) 

The magnetic properties of the conducting fluid are de- 
scribed by the magnetic Prandtl number which is the ratio 
of the kinematic viscosity v to the magnetic diffusivity 77, 
Pm = v/r), /Ltrj denotes the magnetic permeability, and p 
denotes the density. The Reynolds number Re is defined as 
Rc = ^ _1 ili n \/i?i n (i? ou t — i?i n )- The liquid used in the 
PROMISE experiment is the euticetic alloy GalnSn giving 
Pm = 1.4 x 10~ 6 . Thus it is reasonable to solve the MHD 
equations (dimensionless) in the small Prandtl number limit 
(Youd & Barenghi 2006; Roberts 1967; Zikanov & Thess 
1998), 

Bo 



V 2 M 



Ha rot b x 



Bo 



(5) 
(6) 



d t u + (u ■ V)tt = -Vp-f 

V 2 b = - rot(tt x Bo/ Bo 

and div u — 0, div b = 0, where u is the velocity field and 
b is the perturbed magnetic field. 

We simulate the above nonlinear equations for a 2D 
axisymmetric flow in cylindrical coordinates (r,<p,z). For 
details on the numerical method and the boundary condi- 
tions see Youd & Barenghi (2006) and Szklarski & Rudiger 




t [s] u ; [mm/s] 

Fig. 1 Profiles of u z (z, t) at r = R in + 0.6 D forperiodic 
cylinders just above the critical characteristic values: Re = 
1000, Ha = 9.5, j3 = 4. The critical Reynolds number in 
this case is Re c = 842. 



(2006). The cylinders are assumed to be perfectly conduct- 
ing, and the endplates are either conducting or insulating. 
For the latter the pseudo-vacuum approximation is usecQ. 



3 Discussion of the results 

From the point of view of the MRI experiment we are inter- 
ested in obtaining a stable, uniform rotation profile which 
for subcritical characteristic parameters is as close as possi- 
ble to the idealized basic state £!o- On the other hand we 
expect a clear pattern of traveling vortices for supercriti- 
cal conditions. For infinite cylinders with an external axial 
magnetic field and liquids with Pm of our interest the ba- 
sic Couette profile flo is not altered until a critical Reynolds 
number of order O(10 6 -10 7 ) is reached (corresponding to 
a rotation frequency / w 100 Hz). For an instability due to 
the additional toroidal field with (3 = 4 we expect Rc to be 
of order O(10 3 ) (implying / w 0.1 Hz), Ha of order O(10), 
and therefore we search for conditions for which the flow is 
closest to the £!o profile for these parameters (for details on 
the critical values see Hollerbach & Rudiger 2005; Rudiger 
et al. 2005). 

Figure [TJ displays values of the velocity component u z 
measured along the z axis at r = R m + 0.6D for super- 
critical values of rotation and magnetic fields. R[ n — 4 cm, 
^?out = 8 cm, and the physical properties of gallium for 
the viscosity and the magnetic diffusivity are used in order 
to obtain values in physical scales comparable with those 
of the PROMISE experiment. The results in this figure are 
for cylinders with periodic boundary conditions so that the 
profiles are not constrained by end effects and are directly 
comparable with results from linear theory for infinite cylin- 
ders. We notice clear traces of the drifting Taylor vortices. 

1 However, even for copper the assumption for perfect conductors may 
be not very realistic. 
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3.1 Reducing endplate effects 

All undesirable effects induced by the endplates arise as a 
consequence of vertical shears near the boundaries. Thus we 
attempt to reduce the shears by using appropriate boundary 
conditions. Some experiments (e.g. Ji et al. 2004; Noguchi 
et al. 2002) must deal with vast rotation rates since a toroidal 
field is not applied (i.e. f3 = 0), and the rigidly rotat- 
ing boundaries dominate the whole flow (see Hollerbach & 
Fournier 2004). In this case it is necessary to split the end- 
plates into many independently rotating rings (Kageyama 
et al. 2004; Burin et al. 2006). When the rotation rates are 
relatively slow, so that the corresponding Reynolds number 
is of order O(10 3 ), the desired result can be achieved ei- 
ther by allowing the endplates to rotate rigidly and indepen- 
dently (see e.g. Abshagen et al. 2004) or by splitting them 
into two rings which are attached to both cylinders. From 
the technical point of view the latter configuration is easier 
to implement and can be considered as a possible extension 
to the next spiral MRI experiment. 

Firstly we consider a criterion according to which we 
say that the boundary conditions are more suitable. In the 
basic state for subcritical parameters and for the case of pe- 
riodic cylinders, the rotational profile of the fluid is f^o (r) 
and is independent of z, and the magnetic perturbations b 
are zero everywhere. Introducing endplates leads to the de- 
velopment of z and r gradients in velocity, especially close 
to the vertical boundaries where ft(r) from the bulk of the 
fluid must match the imposed conditions at z = 0, z = H. 
Consequently two Ekman vortices, new currents and mag- 
netic fields are generated (we assume the lids to be insulat- 
ing unless explicitly stated otherwise). Any deviation from 
fio will result in generating an azimuthal component of the 
magnetic field, b^, which enters the momentum equation 
(and in our 2D axisymmetric formulation it is the only term 
which gives rise to the Lorentz force). Vertical profiles of 
in the middle, i.e. for r = D/2, D = i? out — R[ n , are shown 
in Fig.|2]for two different boundary conditions. If the end- 
plates rotate rigidly with the outer cylinder, i7 cn d = ^out, 
the Ekman circulation at the bottom lid has a clockwise di- 
rection, if they rotate with the inner cylinder, f2 cnc i = f^m, 
counter-clockwise - all the gradients have opposite sign. We 
conclude that, not surprisingly, there exists a condition with 
flout < ^ondmin < ^in for which the shears are minimized, 
and the generated magnetic field as well. 

We are interested in obtaining a rotational profile for 
which the energy in b^ is minimized, 

E b = J J b4r,z) 2 drdz, (7) 

where the integration is done over the total volume. As a 
measure of the deviation from Hq one could also consider, 
for example, the kinetic energy of the flow, but our aim is to 
obtain good rotational profiles also for Ha of order 10 (and 
(3 = 0). This is not necessarily a good approach since the 
axial field can inhibit the flow velocity while the rotational 
profile will still be significantly different from Slo- 
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Fig. 2 Vertical profiles for induced 6^ in the middle of 

the gap, r = D/2 for (i = 0,Rc = 100, Ha = 1 and 
pseudo-vacuum boundary conditions. f2 en d = ^out> 

fiend — f&in- 
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Fig. 3 The magnetic energy Eb as a function of an- 
gular velocity of independently rotating lid fiend for 
= 0, Rc = 100. Ha = 1, Ha = 10. 

Figure [3] shows how Eb depends on the rotation rates of 
the rigid endplates, Sl nl < £! C nd < ^out- We notice that a 
minimum occurs for f£ en dmin ~ 0.3(f!j n — O ou t) + O ut 
which is even three orders of magnitude smaller than for 

When considering endplates divided into rings, we as- 
sume that the ring attached to the inner cylinder has a width 
w\, the other one attached to the outer cylinder has W2 = 
D — w\. Because it is not obvious which value for w\ should 
be chosen we search for the optimal W\, i.e. for which Ej, 
has a minimum, by performing simulations for several dif- 
ferent values (Wendt 1933, for example, used W\/D = 0.5). 
From Fig. [4] we see that the energy of the induced b^ has a 
minimum for W\JD w 0.4 which is roughly independent 
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Fig. 4 The magnetic energy Eb as a function of radius of 

the inner ring for (3 — 0, Re = 100. Ha = 1, 

Ha = 10. 
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of the applied axial magnetic field. It has also been checked 
that the minimum holds for larger Reynolds numbers (for 
the Fig. [3] as well). Again we notice the improvement of 
Eb of two to three orders of magnitude when compared to 
one end-ring attached either to the inner (wi — D) or outer 
(wi — 0) cylinder. The minimum value is very similar to 
that for independently rotating endplates. 

A qualitative view of the resulting rotational profiles 
gives Fig. [5] which displays deviations of fi(r) - the an- 
gular velocity averaged in the z domain - from f2 (r) for 
different rotational properties of the endplates and for var- 
ied Re and magnetic fields. The case with independently 
rotating endplates refers to boundary conditions where both 
lids rotate with angular velocity fi e nd corresponding to the 
minimum value of Eb- For comparison we also present the 
case for two rings attached to the cylinders with equal width 
Wi = w% = D/2. 

We see that applying independently rotating or split end- 
plates produce significantly more suitable profiles - flatter 
and closer to 1. We also notice that using wi — 0AD gives 
somewhat better results than wi = 0.5D, especially for 
r > R[ n + D/2 where the former profile is almost flat. 

3.2 Influence of the toroidal field 

Figure [6] shows values of the velocity component u z , simi- 
larly like Fig.Q]but for finite cylinders. The velocity field u z 
in the basic state, i.e. = for which there is no instability, 
is u z = everywhere when considering infinite or periodic 
cylinders for Re « 10 3 , Ha * 10. For the enclosed cylin- 
der this is not the case. In Fig. [6^ (left) we present results 
for symmetrically, rigidly rotating (with fi QU t), insulating 
endplates. We notice that u z is quite large and, more im- 
portantly, time dependent (this is even more evident for u z 
closer to the inner cylinder). The right panel in this figure 



a 

la 




2.0 



Fig. 5 Deviations of the averaged Cl(r) from the ba- 
sic state Slo (r) for different vertical boundary conditions; 
rigidly rotating endplates (both with f2 C nd): ^cnd = 

^outs ' ' ' ^cnd — ^iii) — — — r2 cnc [ — ^cndmini divided 

into two rings: ■ - ■ — w-y — 0.5, • • • — • • ■ Wi = 0.4. 
(a) Re = 1000, Ha = 0, (b) Re = 1000, Ha = 10. 

displays the same flow with the toroidal field applied, u z {z) 
is averaged in time and subtracted in order to filter out the 
background. We clearly see the instability and structure of 
traveling vortices, the frequency of this motion agrees with 
the predictions based on the linear analysis (see e.g. Riidiger 
et al. 2005). 



© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Wcinhcim 



www.an-journal.org 



Astron. Nachr. / AN (2007) 



503 



j3 = j3 = 6 




500 1000 1500 2000 500 1000 1500 2000 



t [s] u, [mm/s] t [s] u, [mm/s] 

Fig. 6 The axial velocity u z (z,t) at r = R ln + 0.6D as a function of time t and z. Left: basic state (3 = 0. Right: 
(3 = 6, the averaged u z (t) is subtracted in order to eliminate the background from the velocity field, except in (d). (a): both 
endplates rotate rigidly with f2 e nd = ^out- (t>): both endplates are divided into rings attached to the cylinders; the inner ring 
has the width OAD. (c), (d): the bottom endplate is stationary Slbot = 0, the upper rotates with Q% Q p = ^out- ( a X (b), (c): 
insulating endplates, Rc = 1775 (f2j n = 0.377 Hz), Ha = 9.5. (d): perfectly conducting endplates, Re = 1000, Ha = 10. 
The traveling wave frequency for (a), (b), (c) is respectively f /ft sn = 0.0294, 0.0253, 0.0292 whereas the linear stability 
analysis yields //0; n = 0.0258. 
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As we have shown above, one can obtain a much bet- 
ter basic state for the finite cylinders by dividing endplates 
into two rings. The results for such conditions are presented 
in Fig. |6j). We notice that the background state quickly be- 
comes entirely steady. Naturally the Ekman pumping mech- 
anism is still present in this case, and traces of two Ekman 
vortices can be seen. The flow, however, is laminar. For 
(3 = 6 the pattern of the traveling vortices is clearly more 
regular (cf. Fig.|6p, right). 

When one considers two endplates with different rota- 
tional properties, additional velocity and current gradients 
in the vertical direction arise and disturb further the flow. It 
is clearly seen in Fig. [6J: that disturbances exist in the case 
where the upper endplate rotates with f2 to p = f^out> an d the 
bottom one is fixed, f2bot = 0. The background flow for 
(3 = is highly irregular and time-dependent, especially in 
the middle part of the cylinder, the circulation close to the 
endplates is roughly steady. Nonetheless, the external B$ 
produces, again, a clear periodic motion with a frequency 
corresponding to that of the helical MRI. 

Using conducting boundaries instead of insulating ones 
leads to an increase of the Ekman circulation and the Hart- 
mann current, the latter being drawn from the plates. This 
current is significantly stronger than the current generated 
in the Ekman-Hartmann layer, and therefore we expect that 
an experiment with conducting plates would undergo addi- 
tional problems due to magnetic forces acting on the fluid. 
Let us consider a perfectly conducting endplates with asym- 
metric rotation (again at the top f2 top = ^out an d at the 
bottom Qbot = 0), then there exists an important gradient 
in the radial current which, acting in concert with the axial 
magnetic field, is strong enough to "drag" vortices in the di- 
rection of decreasing field strength. This situation is shown 
in Fig. [6}i where we see a periodic vertical motion even if 
(3 = 0. Moreover, if we introduce a toroidal field with ap- 
propriate sign (i.e. positive in this case) it will act against 
the force due to the current gradient and can reduce the pe- 
riodic vertical motions in the flow (Fig. [6}l, right panel). If 

would have a different sign both effects would interact 
resulting in a highly irregular time-dependent behavior. 

We notice that in the real PROMISE experiment the bot- 
tom endplate rotating with £l out (which, after taking into 
account the directions of rotation and the applied magnetic 
field, corresponds to the top endplate in our simulations) 
was made of copper, and the stationary top endplate (bot- 
tom in the simulations) was made of plexiglass. Therefore 
an additional asymmetry in the magnetic boundary condi- 
tions was present. Although copper is a good conductor it 
should not be directly compared with perfectly conducting 
boundaries used in the simulations since the latter represent 
stronger assumptions and induce stronger currents. How- 
ever it is clear that using insulating material on both ends 
would prevent additional currents from disturbing the flow. 




500 1000 1500 

t [s] u, [mm/s] 

Fig. 7 Profiles of u z (z,t) for Re = 1775, Ha = 9.5, (3 = 
2 and rigidly rotating ends with O e nd = Slout- The criti- 
cal (3 C for the corresponding Re, Ha in the limit of infi- 
nite cylinders is (3 C = 2.56, and one would expect that the 
traveling wave decays. This is not the case for the bound- 
ary conditions shown here where a clear periodic motion is 
visible. Its frequency f/Cl- m = 0.0124 agrees with the pre- 
diction of a linear analysis for marginal stability in the limit 
of infinite cylinders, yielding f/il- m = 0.0120. However, 
the latter approach yields negative growth rate (exact num- 
bers for frequencies and wavenumbers for the linear results 
presented in this paper were provided by R.Hollerbach who 
used them to generate figures in Riidiger et al. (2006). 

3.3 Critical values 

Noting that the background state for sufficiently fast ro- 
tation and rigidly rotating endplates fi e nd = Slout is not 
steady, it is interesting to investigate what happens when a 
spiral magnetic field with strength below the critical value is 
applied. One could expect that a viscid process (like the Ek- 
man pumping) excites fluctuations which could then be am- 
plified and, due to geometry of the applied magnetic field, 
drifting. 

Figure [7] shows that for endplates causing strong dis- 
turbations the traveling wave can indeed be observed even 
for subcritical characteristic values. This is also somewhat 
in agreement with the experiment - traces of moving vor- 
tices were observed for states which are stable in the limit 
of infinite cylinders. We see that the amplitudes of the ver- 
tical component of the velocity u z are almost unchanged 
when compared to the background state (Fig. [6^, left). Al- 
though the pattern of the vortices is not very regular, there 
exists a clear frequency peak for the vertical traveling wave. 
The frequency and the drift direction (which is reversed by 
a sign change of, for example, (3) corresponds to results of 
the linear analysis for infinite cylinders. This leads to the 
conclusion that, although excitations do not grow due to he- 
lical MRI, still the same mechanism is responsible for the 
drift. 

If two rings are used and the basic state is steady the sit- 
uation changes since the additional excitations due to the 
endplates are minimized. Surprisingly, it is possible that 
even for supercritical parameters the traveling wave, al- 
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500 1000 1500 2000 

t [s] u z [mm/s] 

Fig. 8 The velocity u z (z,t) for Rc = 1775, Ha = 
9.5, (3 — 5 and endplates divided into two rings with 
Wi = 0.4. Although in the limit of infinite cylinders the flow 
is unstable, we clearly see that the disturbances (which de- 
veloped after a sudden switch on of the magnetic field) de- 
cay. The frequency of the decaying wave is / /Qi n = 0.0237 
which again is in agreement with / from the result of the 
linear analysis, f/Q- m — 0.0232. 

though excited for a moment, decays (see Fig. [8}. It is still 
possible to get sustained instability by increasing, for exam- 
ple, (3 (see Fig.|6p). 

The reason for this damping might be height of the 
cylinders which does not match an integer value of the ver- 
tical wavenumber k. For Re = 1775, (3 = 5, Ha = 9.5 the 
corresponding wavelength is A = 27r_D/2.1728 and does 
not suit the assumed aspect ratio T = H/D = 10. If the 
height was changed to T = 4 A = 11.57 the observed decay 
of the traveling wave was significantly slower, so slow that 
after the sudden switching on of the external azimuthal mag- 
netic field the wave could be observed with the PROMISE 
facility still several hours later. Bearing in mind that wave- 
lengths for given Reynolds numbers are longer with de- 
creasing beta (for = 3, k = 1.45-D -1 , = 1, A; = 
0.6Z3 -1 ), the constant height of the cylinders (T = 10 in 
the experiment) can be an issue when looking for critical 
numbers. It should be noted that due to the boundary lay- 
ers the effective region where the traveling wave can exists 
for configuration with two rings is smaller than T by ap- 
proximately 0.5D (distances up to about Q.25D from the 
endplates are influenced by their presence). 

Although endplates clearly can serve as the source of 
viscous excitations and the axially traveling wave develops 
also for subcritical parameters, we shall notice that there are 
no periodic motions in the background state. In this sense 
the "imperfect" background state serves as a catalyst for the 
helical MRI instability. When the endplates are divided into 
rings the resulting hydrodynamic flow is laminar, and only 
after the magnetic field is applied the periodic fluctuations 
occur, and, moreover, their frequency corresponds exactly 
to that predicted from the linear analysis for infinite cylin- 
ders. 




500 1000 1500 

t [s] u, [mm/s] 

Fig. 9 The velocity u z (z,t) for Rc = 1775, Ha = 
9.5, (3 = 6 and endplates rotating differentially so that 

"end = ^0- 

Liu et al. (2006) suggested that the observed fluctuations 
can have theirs origin in the underlying hydrodynamical un- 
steady flow as reported, for example, by Kageyama et al. 
(2004). In the latter work the purely hydrodynamic flow for 
Re w 1000 with rigidly rotating ends fi cn d = ^out and 
short aspect ratio r = 1 was already unsteady. We confirm 
these results with the method used here. However if longer 
cylinders are used, like T = 10, the flow becomes steady 
for Re = 1000, and only after imposing strong enough mag- 
netic fields (say Ha = 12, (3 = 6) a traveling wave develops 
with a frequency that matches calculations from the linear 
analysis. 



3.4 Differentially rotating ends 

We have also performed simulations for differentially ro- 
tating plates with ideal Couette profile for periodic cylin- 
ders so that f2 cn d (r) = ^oM- I n another recent work (Liu, 
Goodman & Ji 2007) it has been shown that for parameters 
corresponding to Rc = 1775, Ha rj 10, (3 ss 4 the travel- 
ing wave decays for such boundary conditions. We confirm 
this result with our method, although our treatment of the 
magnetic boundaries is simplified. 

The explanation for this fact might be again the inappro- 
priate height of the cylinders which is far from an integer 
value of the expected vertical wavelength. For these param- 
eters A = 3.476D according to the linear theory so that less 
than three wavelengths can fit in the container. On the other 
hand, if (3 = 6 is used, A = 2.49L» and then r = 10 al- 
most exactly corresponds to 4A. From Fig.|9]we see indeed 
that in this case persistent fluctuations exist with a frequency 
corresponding to the helical MRI instability. We have also 
made calculations for (3 = 4 with longer cylinders so that 
H = 4XD = 13.90£> and T = 5XD = 17.38D. In each 
case a sustained traveling wave has been observed. It should 
be mentioned that the vortices do not develop very close to 
the upper boundary so that it is convenient to take longer 
cylinders. 
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4 Summary 

It is easier to perform an experiment showing spiral MRI 
because a much slower rotation of the cylinders is required 
for the instability to set in compared with MRI with an axial 
magnetic field only. Moreover, there exists additional quan- 
tity, i.e. drift frequency, which is easy to measure and can 
serve as an important indicator for the associated phenom- 
ena. It is claimed that in the PROMISE experiment frequen- 
cies and amplitudes corresponding to the spiral MRI were 
observed, and the results agreed with theoretical calcula- 
tions of both linear and nonlinear 2D simulations (see Ste- 
fani 2007) for a review). However it is still possible to im- 
prove the experiment so that the basic state is a completely 
steady flow. 

In this paper we have presented a relatively simple and 
inexpensive modification which is suitable for such an im- 
provement. Firstly, the endplates should be both made of 
insulating material and both should rotate in the same way 
so that the system is symmetric in the z direction. Secondly, 
it is convenient to divide the lids into two rings which can 
be attached to the cylinders so that no separate driving is 
needed. The optimal, width of the inner ring, in the sense of 
minimizing the induced azimuthal magnetic field, is 1.6 cm 
for the current experimental setup. 

Our calculations also show that spiral MRI modes can 
be driven by endplate effects even for subcritical charac- 
teristic values (see Fig. 0. On the other side when provid- 
ing a steady background flow by applying rings one has to 
pay more attention to the height of the cylinders and to take 
into account the vertical wavelengths of the traveling wave 
which depend on the magnetic configuration. For the cur- 
rent aspect ratio r = 10 and Re = 1775 it is reasonable 
to consider Ha = 9.5, = 6 which almost exactly corre- 
sponds to r = 4A. 
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